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Abstract 

The dynamics of three mutually coupled cortical neurons with time delays in the coupling are 
explored numerically and analytically. The neurons are coupled in a line, with the middle neuron 
sending a somewhat stronger projection to the outer neurons than the feedback it receives, to 
model for instance the relay of a signal from primary to higher cortical areas. For a given coupling 
architecture, the delays introduce correlations in the time series at the time-scale of the delay. It 
was found that the middle neuron leads the outer ones by the delay time, while the outer neurons 
are synchronized with zero lag times. Synchronization is found to be highly dependent on the 
synaptic time constant, with faster synapses increasing both the degree of synchronization and the 
firing rate. Analysis shows that pre-synaptic input during the inter-spike interval stabilizes the 
synchronous state, even for arbitrarily weak coupling, and independent of the initial phase. The 
finding may be of significance to synchronization of large groups of cells in the cortex that are 
spatially distanced from each other. 

PACS numbers: 
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I. INTRODUCTION 



There is a significant amount of research showing that spike coincidence of neurons en- 
codes information. Synchronization of neural activity has been involved in such important 
phenomena as ability to recognize objects (by binding different attributes) [9|, olfactory dis- 
crimination 2], and has even been proposed as one of the neural correlates of consciousness 

n 

[3|. Abnormal neural synchronization has been linked to the onset of epileptic attacks. 

Evidence of long range correlation leading to synchronized clusters has been observed in 
the visual cortex [5j where synchronized oscillatory responses are conjectured to contribute 
to the binding of distributed features in a visual scene. Further research in on the cat 
visual cortex supports that oscillatory activity and firing patterns can contribute to long 
range synchrony among cortical cluster of neurons. Other mechanisms involving gamma 
oscillations are also conjectured as one underlying cause for long range correlations in the 
human brain. 

On the modeling front, recent models involving self-feedback loops in simple motif net- 
works have shown that coupling architecture may improve long range correlations. Traub 
et al [3] uses a linear self-feedback loop in the coupling architecture of cortical neurons 



to inhance synchrony. In [10|], theory and experiments on semi-conductor lasers, which 
exhibit bursting phenomena similar to neurons, show that complete synchronization may 
occur with self-feedback loops. In addition, lasers coupled in a line without self-feedback 
have been shown to synchronize the end lasers Q] , leading to the possibility of long range 
correlations with delay. 

Synchronization of nearby cells is often the result of receiving common input, such as 
when a pyramidal cell sends projections to a targeted area in the cortex. While pyramidal 
cells tend to target specific areas, matrix projection cells from the thalamus reach in a diffuse 
manner into adjacent cortical areas helping to synchronize the activity of large populations 
of cells. Since many strong synaptic connections between different brain areas tend to be 
one-way, many of the networks in the brain possess a hierarchical structure, which to a first 
approximation is a forward network, modulated by a weaker modulatory feedback ([4f)- For 
example, there is a strong forward projection from the thalamic LGN area to VI or from VI 
to MT, with weaker modulatory feedback projections that modulate the magnitude of the 
cell's response (|14|). It has been suggested that exessively strong mutually coupled loops 
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in the brain would promote uncontrollable oscillations, such as in epilepsy (|13|). 

With this in mind, we propose a model of synchronization of cortical cells receiving a 
common time-delayed input from a different area of the brain, such as the case with projec- 
tions of pyramidal cells to other cortical areas or the thalamus to the cortex, and receiving 
a weaker time-delayed feedback. The central dynamical question addressed considers when 
the synchronous behavior of such a neural network is stable, particularly in regard to the 
synaptic coupling strengths, and the synaptic time constant. We find that shorter synaptic 
time constant of the target cells promotes synchronization, and even increases the firing 
rate, for the same strength of input. 

The paper is organized as follows: In Section II, the basic model is set up, (see Fig. [Q 
and presented along with numerical results that show the dependence of synchronization 
and firing rate on the strength of coupling and synaptic conductance. Section III, analyzes 
the synchronization observed for short synaptic times by linearizing about the dynamics of 
two nearby trajectories. It shows that fast synaptic input tends to synchronize two nearby 
trajectories, regardless of their phase with respect to the input signal. Section III concludes 
and summarizes. 



II. BASIC MODEL AND NUMERICS 

Neurons can be largely divided into two classes, dependent on their spiking properties. 
Class I neurons can be stimulated to fire at an arbitrarily low frequency, due to a saddle-node 
bifurcation, with increasing frequency as the magnitude of the stimulating current increases. 
Class II neurons, on the other hand, only begin to fire at relatively high frequency, with 
their limit cycle resulting from a subcritical Hopf bifurcation. Class II neurons are well 
represented by the squid axon, while a large majority of the mammalian neurons are of 



the Class I type. Dynamics of a human neo-cortical neuron in the absence o 
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synaptic 



connections are well approximated by the following Class I neuron equations: 

^ = -{17.81 + 47.58V + 33.8V 2 } (V - 0.48) - 2QR (V + 0.95) + I = F (V, R) + I 
dH 1 

— = — (-R + 1.29V + 0.79 + 3.3 (V + 0.38) 2 ) = G (V, R) , (1) 
dt t r 

where V and R are voltage and recovery variables, respectively, and tr = 5.6 ms. The dR/dt 
equation is written as a sum of a linear term, for the normal Na + and K + currents, and 



a quadratic terms in V to approximate the the transient potassium current contributions 



171 ] . The above model has been optimized to provide an accurate quantitative fit to the 



shape of a regular spiking neuron potentials obtained from human neocortical neurons [19 ]. 
Figure [2] shows the limit cycle of a cortical neuron in Eq. ([T]) when the applied current, /, 
is above the bifurcation value, resulting in a saddle-node bifurcation. Due to a quadratic 
term in the recovery variable, the spike rate of a cortical neuron can be arbitrarily low, 
for low currents, and increases as I is increased. The above neurons can be coupled by 
adding an additional term to dV/dt in Eq. (pQ) proportional to g (V — E syn ), where g is the 
synaptic conductance variable, E syn = for excitatory synapses and —0.92 for inhibitory. 
The synaptic conductance, g, is obtained from the following equations, commonly used for 
synaptic coupling [l2], 

-77 = ( — / + Hstep (Vpre ~ ^)) 

Tsyn 

| = ^-(-» + /) P) 

UiL 1 syn 

where H step (x) — 1, if x > and zero if x < 0, V pre is the voltage of the presynaptic neuron 
and r syn is the synaptic conductance. In numerical simulation, Q = —0.20 mV was chosen 
121 ] . The reason for using two synaptic equations is that, depending on r syn , the conductance 
will peak after V pre , continuing to depolarize the membrane after the end of the presynaptic 
spike (see Figure [3]). This type of response is consistent with physiological data. For brief 
stimulus spike at t = 0, the two synaptic equations produce a response proportional to [kJ 

9 = (t/ T 2 syn ) exp (-*/ T syn ) (3) 

Figure [3] shows g for different values of r syn . In each case, the area under the curve is the 
same (equal to one), with r syn determining the width of the post-synaptic spike. 

Using Eqs. (JT]) and (j5J) we can now set up the basic model of three mutually coupled 
neurons. The basic set-up is shown in Figure [1] where 3 neurons are coupled in a line. The 
middle neuron sends an excitatory time-delayed signal to the two outer neurons which are 
coupled to the middle one via weaker delayed coupling with a longer synaptic time constant. 
This model reflects the typically observed hierarchy in the brain (described in the introduc- 
tion), where the neurons in one brain area, such as the thalamus, send strong excitatory 
connections to cortical neurons, receiving weak modulatory feedback. Since the equations 
of the outer neurons are identical, receiving the same synaptic input from the middle cell, 
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their synchronization would indicate that a group of identical neurons fire synchronously 
when subjected to a particular synaptic input. Using Eq. ([I]), the equations for the coupled 
neurons are given by: 

^ = F (Vi, Ri) 5 i9i (V, - E syn ) 

^ = (4) 

where % is the index of each neuron, with outer neurons receiving the same current, = I s , 
and the same synaptic strength, 5\ = 5^ = 5. The current input to the middle neuron 
is higher than to the outer ones: I2 > I s , leading to the lower uncoupled spiking rate for 
the outer cells. This was done so that the higher activity of the inner cell drives the outer 
ones, as might be the case in a typical hierchical network. The synaptic input to the middle 
neuron is weaker, 5 2 < 8. The synaptic time-constant in Eq. ([2]) is the same for the outer 
neurons, r syn \ = r syn3 = T syn , but longer for the middle neuron, to model the affect of a 
slowly- varying modulatory feedback. The conductances, gi, are obtained from Eq. (T5]), with 
presynaptic voltage to the middle cell given by: V pre 2 = V\{t — r d ) + V3 (t — r d ), and to the 
outer cells by V pre \^ = V 2 (t — r d ). The delay in the propagation of the signal is given by r d . 



The effect of delays on correlation and synchronization 

While the length of the synaptic delay, r d , does not seem to affect the degree of synchro- 
nization, it has a substantial affect on correlations and phase relations between neurons. 
Figure H] shows correlations in spiking output between the middle and outer neuron, and 
between outer neurons for r d = 10. The x-axis indicates the time-shift at which the correla- 
tions function was computed. The outer neurons are synchronized, since C\z = 1 at t = 0. It 
can be seen that the greatest correlations between the inner and outer neurons occur when 
t is shifted by the delay time, r d . The delay also creates spikes in correlation at intervals 
of twice the delay time, as can be seen in Fig. HI Since, in this case, the outer neurons are 
synchronized, the spikes indicate that there are self-correlations in the time series of a single 
neuron at intervals of 2r d . This is the round-trip, or feed-back, time, since its the minimum 
time that it would take for a signal to travel from one of the neurons, affect the target, 
and get back to that same neuron. It follows that time delays lead to self-correlations in 
the spike trains that are not observed when delays are absent, and which may lead to more 
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regular patterns in the time-series data. 

For reasons explained in the introduction, 5 > 62 and r syn < r syn 2 were used to model the 
affect of a stronger forward and a weaker modulatory feedback. For this type of coupling, 
the middle neuron leads the outer ones by the delay time, r^. This type of phase-locking 
behavior has been observed in other time-delay systems, such as lasers 
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2 OH. Thus the 



time-shift between two correlated spike trains should be directly related to the delays in 
transmission, with the input leading the output by the delay time. 



The effect of the coupling strength on firing rate and synchronization 

Increasing the synaptic coupling in general increases the firing rate of the neurons for 
excitatory synapses. Above a certain value of the coupling strength, there is phase-locking 
between the inner and the outer neurons, whereby all firing rates are equal. The bifurcation 
value of synaptic coupling at which 1 : 1 frequency locking occurs depends on the current 
injected into each cell. In the absence of coupling, a significant difference in injected current 
for the outer and inner neurons, (J s = 0.22 and I2 = 0.5 in simulations) leads to a big 
difference in firing rate. Figure [5] shows a typical voltage trace for the uncoupled, 5 = 0, 
and coupled neurons, 5 = 4, where coupling is sufficiently strong to cause phase-locking. 
Phase-locking or 1 : 1 frequency locking occurs for 5 > 3.4. This type of behavior has 



been observed for mutually coupled neurons in the absence of delays 12j. For lesser values 
of the coupling strength, different frequency locked behaviors are observed, with frequency 
ratio between 1 and fi/ f , where fa and /„ are the frequencies of inner and outer neurons, 
respectively, in the absence of synaptic coupling. Figure [6] shows 1 : 4 frequency locking 
that occurs for weak synaptic coupling: 5 = 0.5, 52 = 0.3. When the coupling is too weak, 
the outer neurons become desynchronized, as shown in Figure The bifurcation value of 5, 
however depends on the synaptic time constant, r syn , of the outer-neurons. Thus Figs. Eland 
[HI which have a very short time constant of r syn = 0.03, shows synchronization at a lower 
coupling strength of 5 = 0.5, compared to a minimum 5 = 1.03 needed for synchronization of 
neurons in Fig. [71 where r syn = 0.5 (a more realistic value for fast synapses). This sensitive 
dependence of synchronization on the coupling strength and synaptic time constant of the 
outer-neurons is explored analytically in Section III. 
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The effect of synaptic time- const ant on synchronization and firing rate 

Figure [9] shows an increase in firing rate of the outer neurons as their synaptic time 
constant, r syn , is decreased. This increase in firing rate may be surprising, since the time 
constant only controls the width of the conductance spike and not the area under the curve, 
as shown in Fig. [31 Thus the contribution of a pre-synaptic spike to a change in post- 
synaptic voltage during an inter-spike interval is largely independent of the synaptic time 
constant. This can be seen by integrating the g (Vi — E syn ) term in Eq. over the interval 
of conductance change. Figure [TOl shows a fluctuation in conductance for a system given by 
Eqs. ([2]) and (j3J), with r syn = 1 and T syn2 = 3. A lower firing rate that occurs for slower 
synapses may be the result of the decay of any increase in voltage during the inter-spike 
interval back to the limit cycle trajectory (see the bottom of Fig. [8]). A pre-synaptic spike 
from the inner neuron can trigger a spike from the outer one, when it is delivered toward the 
end of the inter-spike interval. Thus a more narrow jump in conductance and the resultant 
jump in voltage, Vi^, may mean that there is less decay before the critical threshold is 
reached, thereby increasing the likelihood of a spike. 

The most noticeable affect of the decrease in r is greater synchronization. This is shown 
in Fig. [TTJ where the outer neurons become progressively synchronized as r syn is decreased 
from 0.5 to 0.2. It can be seen that for relatively weak coupling of 8 — 1, 82 — 0, and 
Tsyn — 0.2, the outer neurons are completely synchronized, after the transients die out. The 
next section analyzes the affect of a presynaptic spike on a fast synapse (small T syn ) in 
synchronizing two nearby trajectories of the outer neurons. 

III. ANALYSIS OF SYNCHRONIZATION FOR WEAKLY COUPLED, FAST 
SYNAPSES 

The three neuron model described by Eqs. (j2J) and (jlj) possesses internal symmetry, e.g. 
its equations of motion do not change if the variables {Vi, Ri} and {V3, -R3} are interchanged. 
It follows that the synchronized regime, where the symmetric variables are exactly equal: 
{V\ — V3 — Vo, Ri = R3 = Rq}, is a solution For an uncoupled system, this solution 
would not be a stable one, since any perturbation along the limit cycle would result in a 
phase-difference. A spiking input, V2, from the middle neuron stabilized the synchronous 
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state when the synaptic time- const ant, r syn is sufficiently short. To study the affect of a 
single pre-synaptic spike on two identical neurons, with nearby trajectories along the limit 
cycle, introduce new variables: V = V% — V3 and R = R\ — R3. These new variables 
correspond to a perturbation transverse to the synchronized state: {Vo,i?o}- Using Eqs. 
([T]), ([3]) and fll]), the linearized dynamics are: 

^ = -n(V )R-h (Vo, R )V-S (t,lT% n ) exp (-i,/r s „„) V (5) 

where n (Vo) = 26(V + 0.95), h(V ,R ) = (-5.03 + 101.4V 2 - 32.45 V + 26i? ), and 
m (Vo, Rq) = — (3.8 + 6.6V0). t s denotes the time since the last spike of the center neu- 
ron. Equation (JS]) is valid for a short synaptic time constant, when only the effect of the 
last presynaptic spike is significant (see Fig. [3]). During an inter-spike interval of the outer 
cells, n, h and m are all positive. This can be easily confirmed by using Fig. [21 where 
—0.8 < V < —0.6 and .18 < R < 0.4 during the inter-spike interval, and calculating the 
lowest possible values of n, h and m for a given range of Vo and Rq- Since {V^, R} denote the 
difference between the two nearby trajectories along the limit cycle, there is a relationship 
between the two variables given by 

R(t) = -l(V )V(t) (7) 

where / is the negative of the slope of the limit cycle at {V (t), R (t)}. From Fig. [21 /(Vo) ~ 2 
during the inter-spike interval. 

A change in V following a post-synaptic spike can now be calculated by substituting Eq. 
© into Eq. ([5]), dividing by V and integrating over At, the duration of the conductance 
spike, (t s /Tg yn j exp (—t s /T syn ). For a short synaptic time constant, r syn < 1, the width of 
At is small compared to the time-scale of neuronal dynamics during the inter-spike interval 
of the outer cells. A change in V following a post-synaptic spike is now given by, 

V(t + At) = V u (t + At) ■ e~ 5 (8) 

where V u (t + At) is the voltage difference of two nearby trajectories at time t + At in the 
absence of synaptic input: V u (t + At) m V (t) exp (At {nl — h)), with n and h defined after 
Eq. §6§ and I given in Eq. ([7]). For small At, so that | At (nl — h) \ <C {1, 6}, most of the 
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change in voltage during a post-synaptic spike is caused by the spike itself so that Eq. (JHJ) 
can be approximated as 

V (t + At) « V (t) e- 5 (9) 

Using Eqs. ©-([ED, the dynamics of R immediately following a narrow post-synaptic spike 
can be approximated as 

dR _ If,, fm(V ,R )\ s 



=f = — -1 + u; ; w e-M R (10) 

dt 5.6 V \ I (V ) J J K J 

From Eqs. (jSj) and a single presynaptic spike from a center neuron acts to decrease the 
voltage difference, V, of the outer neurons by a factor of exp (—5). This presynaptic spike 
also decreases R by decreasing the positive contribution from the m/l > term in Eq. (flOl) . 
Thus the effect of the presynaptic input is to decrease the perturbation from synchronized 
state of the outer neurons, pushing their trajectories closer in phase-space. 

Eqs. (l8"j)- ([TUj) show that a spiking input from the middle neuron has a stabilizing affect 
on the synchronized state when the synaptic time constant is short. Since the difference 
in trajectories, V, R is taken along a limit cycle, the maximum Lyapunov exponent in the 
absence of synaptic coupling would be zero (corresponding to the displacement along a 
trajectory in phase-space), and the transverse exponents must be negative since the systems 
is dissipative. It follows that, for sufficiently small T syn , a common synaptic input acting 
during the inter-spike interval should eventually synchronize the neurons, even for weak 
synaptic coupling. 

Figure [8] shows the synchronization of the outer neurons for a very short synaptic time 
constant, r syn = 0.03 and weak coupling, 5 = 0.5. After the transients die out, the outer 
neurons become synchronized, thereby falling on a straight line in the V\ vs V3 plot. As 
can be seen in Figure El for weak coupling, there is a big difference in firing rate between 
the outer and the middle neuron, due to differences in injected current. It follows that the 
synchronization is not due to phase-locking between the middle and the outer cells. The 
effect of the spiking input on the limit cycle trajectory can be seen at the bottom of Figure 
[8]) . There is an integer ratio between the inner and the outer frequencies, whereby the onset 
of a spike in the outer neuron is triggered by spiking input from the middle one, delivered 
toward the end of the inter- spike interval. The phenomena is similar to the subharmonic 
resonance where the limit cycle responds at a subharmonic of the stimulus frequency 12]. 
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Eqs. (IH|)- (|TO|) show a strong dependence of synchronization on the coupling strength. 
This is confirmed by numerics. Figure [7] shows V\ vs V3 for two slightly different values of 
5,5 = 1.02 and 5 = 1.03. A slight change in 5 leads to an onset of synchronization between 
the outer neurons. 

IV. CONCLUSION 

A model of three mutually coupled cortical neurons with delays was studied using analysis 
and numerical simulation. The outer neurons were stimulated with smaller current and had 
a much lower frequency in the uncoupled case, with their frequency significantly increasing 
depending on the strength of synaptic coupling with the middle neuron. At higher values 
of the synaptic coupling constant, typical phase- locked behavior and 1 : 1 frequency locking 
was found between the middle and the outer neurons, with different frequency locking ratios 
as the synaptic strength was lowered. It was found that delays affected the time-series 
data by introducing correlations at the time-scale of the delay. While the spiking behavior 
in the synchronized case was fairly regular, this effect would be interesting to explore for a 
more complicated, chaotic spike train that can be achieved by incorporating slow adaptation 
currents into the neuron model. In the case of phase or frequency locking, the middle neuron 
leads the outer by the delay time, 

While synchronization of outer neurons was sensitive to the synaptic strength, the synap- 
tic time constant for outer neurons, r syn was also highly significant. It was found that shorter 
synaptic constant substantially improves synchronization, able to synchronize the neurons 
even when the coupling strength is very weak. A short synaptic constant was also able to 
induce 1 : 1 frequency locking at a much weaker mutual coupling than this type of dynam- 
ics might be expected, for the input currents used. Analysis of dynamics for fast synapses 
showed that spiking synaptic input during the inter-spike interval stabilized the synchroniza- 
tion manifold, even for arbitrarily weak coupling, and independent of the phase relationship 
between the inner and outer cells. This indicates that even a very weak synaptic input can 
synchronize cells, as long as the synaptic time constant is sufficiently short. The finding 
may have significance in synchronizing large groups of cells in the cortex via weak synaptic 
input from other areas, such as the thalamus, or other areas in the cortex proper. 



10 



ACKNOWLEDGMENTS 



This work was supported by a grant from the Office of Naval Research. ASL is currently 
a post doctoral fellow with the National Research Council. 



Electronic address: jalandsma® cantor. nrl. navy. mil^ 



[1] A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization: A universal concept in nonlinear 

science (Cambridge University Press, Cambridge, 2001). 
[2] M. Stopfer, S. Bhagavan, B. H. Smith, and G. Laurent, Nature 390, 70 (1999). 
[3] F. C. Crick and C. Koch, Sci.Am. 273, 84 (1995). 
[4] F. C. Crick and C. Koch, Nature 391, 245 (1998). 

[5] A. K. Engel, P. Konig, A. K. Kreiter, and W. Singer, Science 252, 1177 (1991). 

[6] I. Fischer, R. Vicente, J. Buldu, M. Peil, C. Mirasso, M. Torrent, and J. Garcia-Ojalvo, 
http: / /arxiv.org/abs /nlin /0605036 (2006) . 

[7] A. S Landsman, I. B. Schwartz, http://arxiv.org/abs/nlin/0609047 (2006). 

[8] J. Guckenheimer, and P. Holmes, Nonlinear Oscillations, Dynamical Systems and Bifurcations 
of Vector Fields (Springer- Verlag, New York, 1983). 

[9] W. Gerstner, and W. Kistler, Spiking Neuron Models (Cambridge University Press, Cam- 
bridge, 2002). 

[10] E. Klein, N. Gross, M. Rosenbluth, W. Kinzel, L. Khaykovich, and I. Cantor, arXivxond- 
mat/0511648. 

[11] P. Konig, A. K. Engel, and W. Singer, Proceedings Of The National Academy Of Sciences Of 
The United States Of America 92, 290 (1995). 

[12] H. Wilson, Spikes decisions and actions (Oxford University Press, New York, 1999). 

[13] C. Koch, Biophysics of Computation (Oxford University Press, New York, 1999). 

[14] C. Koch, The Quest for Consciousness (Roberts and Company Publishers, Englewood, Col- 
orado, 2004). 

[15] E. Rodriguez, N. George, J. P. Lachaux, J. Martinerie, B. Renault, and F. J. Varela, Nature 
397, 430 (1999). 

[16] R. D. Traub, M. A. Whittington, I. M. Stanford, and J. G. R. Jefferys, Nature 383, 621 



11 



(1996). 

[17] R. M. Rose, J. L. Hindmarsh, Proc. R. Soc. Lond. B 237, 267 (1989). 

[18] T. Heil, I. Fischer, W. Elsasser, J. Mulet, and C. R. Mirasso, Physical Review Letters 86, 795 
(2001). 

[19] R. CFoehring, N. MLorenzon, P. Herron,and C. J. Wilson, J. Neurophysiol. 66, 1825 (1991). 
[20] J. K. White, M. Matus, and J. V. Moloney, Physical Review E 65 (2002). 



12 




FIG. 1: The basic model of three synaptically coupled neurons. The strength of the coupling from 
the middle to the outer neurons is given by 5 and from the outer to the middle by 62 ■ The injected 
currents are shown. All synaptic coupling has propagation delay of and a synaptic time constant 
of T S yn a n d T S y n 2 for the outer and the middle neurons, respectively. 
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FIG. 2: Limit cycle of an uncoupled cortical neuron, given by Eq. (pQ). / = 0.5 nA. 
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FIG. 4: Correlations of spiking output. Top: between the outer and the middle neuron. Bottom: 
Between outer neurons, = 10 ms, 5 = 4, 82 = 2, T syn = 1 ms, T syn 2 = 2 ms, = I s = 0.22 nA, 
I 2 = 0.5 nA. 
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FIG. 5: Left: spiking in the absence of coupling. Right: coupling, 5 = 4, 62 = 2, r syn = 1, 
T S yn2 = 2. In both cases, I s = 0.22, I2 = 0.5. 
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FIG. 6: Spiking voltage for weak coupling, 5 = 0.5, r syn = 0.03, 82 = 0.3. There is 1 : 4 frequency 
locking between the outer and the middle neurons. 
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FIG. 8: Synchronization for weakly coupled, fast synapses. 5 = 0.5, T syn = 0.03, 82 = 0.3. Top: 
synchronization after the transients die out. Bottom: Limit cycle of one of the outer neurons. 
Same parameters as in Fig. [6l 
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FIG. 9: Increase in firing rate as T syn decreases. Left: T syn = 2. Right: T syn = 0.1. In both cases, 
5 = 1, 5 2 = 0.2, T syn2 = A,I S = 0.22, J 2 = 0.5. 
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FIG. 10: Conductance for T syn = 1, T syn 2 = 3. All else as in Fig. H] 
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